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Abstract. We use a simple gas model to study non- 
equilibrium aspects of the multiparticle dynamics rele- 
vant to heavy ion collisions. By performing numerical 
simulations for various initial conditions we identify sev- 
eral characteristic features of the fast dynamics occurring 
in implosion-explosion like processes. 

Introduction. Thermodynamic equilibrium at an in- 
termediate stage of a heavy ion collision is frequently a 
basic assumption in models of the colliding nuclear mat- 
ter. These models range from statistical models of nu- 
clear multifragmentation to the fluid dynamical models 
of the quark gluon plasma. In contrast, microscopic mod- 
els of molecular dynamics type (e.g. RQMD, FMD and 
NMD), which are based upon constituent interactions do 
not have this assumption built in. Such models are ap- 
propriate for testing to what extent thermodynamic equi- 
librium is actually achieved. 

Previously we have used the Nuclear Molecular Dy- 
namics (NMD) model [Q to describe the collision dy- 
namics and clusterization process in nucleonic matter. 
In a special study of the time evolution of a finite sys- 
tem of colliding particles we noticed that the phase space 
distance d between particles in the combined coordinate 
and momentum space has a distinct time evolution for 
different groups of particles: For nucleons which form 
clusters, d is initially large and rapidly decreasing until it 
stays small and rather constant from a certain time. This 
can be interpreted as freeze-out of the hard core scatter- 
ing interaction. This interpretation was supported by the 
fact that for interactions with no hard core scattering the 
phase space distance for clusters always remained small. 
Averaged over all nucleons the phase space distance d 
grows with time reflecting the overall expansion of the 
system. 

The aim of the present paper is to study the behavior 
of a finite piece of exploding matter before and around 
freeze-out of hard core interaction. Before and at this 
stage, the evolution of the matter is highly dynami- 
cal, leading to characteristic space-time distributions of 
freeze-out and flow. Whether the dynamics actually re- 
sults in thermodynamical equilibrium before frcczc-out is 
not clear a priori. Different aspects of this problem have 
been studied previously within various approaches [^]-|| . 
In a previous paper Q we employed a new variable, the 
pseudo entropy, to study the equilibration. In the present 



letter we use theoretically well-founded concepts to inves- 
tigate the dynamics of the exploding matter. 

Classical gas model. To this end we study the dy- 
namical properties of a simple gas model, where the con- 
stituent particles interact like billard balls by classical 
Newtonian dynamics. A first application of this model 
of nuclear collisions was made in . We consider a gas of 
identical classical balls of radius r c . They perform classi- 
cal non-relativistic hard-sphere elastic scatterings at the 
impact parameter b < 2r c with conservation of energy, 
momentum and angular momentum. (The balls do not 
have intrinsic rotation.) The initial configuration consists 
of A such particles placed randomly within a sphere of 
radius R, rejecting configurations where particles overlap 
within the hard core distance. This puts a constraint on 
R which should not be too small. 

We prepare the system at certain energy and density, 
and let it evolve until the interactions between the con- 
stituents cease. We further impose spherical symmetry 
on the system. In this way we minimize specific compli- 
cations of the interactions and the geometry, and we can 
concentrate on a few basic properties of the behavior of 
small shortlived systems, which we believe are of general 
interest for heavy ion physics. 

We study the evolution of the matter in the gas ball in 
a simple dynamical process, which is initiated by ingoing 
flow of the particles (implosion). In this situation the sys- 
tem will first compress and then expand, thus simulating 
some basic features of a heavy ion collision. 

The initial velocities of the particles are chosen as a 
superposition of thermal and collective motion. We use 
a spherically symmetric linear profile for the initial col- 
lective velocity field 



v{r) = v Qf r/R 



(1) 



where Vof is a model parameter. In our simulations we 
fix the total energy E = Efi + Eth, and vary the fraction 
77 of the flow energy, 77 = Efi/E, where Efi and E t h 
are the flow energy and the thermal energy, respectively. 
Because of the way in which the system is built up, these 
energies will fluctuate from event to event with a relative 
uncertainty of order A~2 . 

We use conventional nuclear scales: the mass of the 
constituent particle is 940 MeV, the hard core radius is 
r c = 0.4 — 0.5 fm, and the initial radius of the gas sphere 
is R = 1.2 A 1 / 3 fm. We have performed simulations with 
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A = 50 and A = 100 varying the flow velocity parameter 
vof from —0.5 to 0.5 (this could be in units of the velocity 
of light, c = 1, but the velocity scale is in fact arbitrary 
in this model), and varying r\ from (thermal explosion) 
to 1 (implosion). 

Freeze-out. A very important characteristics of evolv- 
ing matter is the scattering rate v{t), defined as the total 
number of scatterings per particle per unit time. It is 
obvious that in an implosion v(t) first grows and then 
decreases, while in a thermal explosion it decreases. In 
Fig. [I] shows v(t) for the two cases. The scattering rate 
roughly follows the matter density, because the mean 
free path A is inversely proportional to the density n 
(A = 1/na, a being the scattering cross section). 
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FIG. 1. The total scattering rate for a thermal explo- 
sion (solid curve) and an implosion (dashed curve) with 
«o/ — —0.5. Here, and in Figs. |2| and ^ r c = 0.5 fm. 

One of the most fundamental questions in heavy ion 
physics is when can a collection of particles be considered 
as coherent matter? In our simplified model we shall con- 
sider particles which collide with other particles as the 
matter, whereas particles which do not scatter any more 
are decoupled from the matter, or frozen out. Therefore 
we study the space-time distribution of the last scatter- 
ing points (the freeze-out field), because this defines the 
space-time limit of the coherent matter. 

Two examples of a freeze-out field are shown in Fig. [|. 
It is seen that the distribution is rather spread, in partic- 
ular for the thermal explosion. This disagrees with the 
Cooper-Frye picture pL which assumes a sharp freeze- 
out hypersurface. A similar conclusion was reached in 0] 
for microscopic simulations of relativistic heavy ion col- 
lisions. 
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FIG. 2. Freeze-out fields for a thermal explosion (w = 0, 
top panel), and for an implosion (77 = 1, lower panel), with 
100 particles. Vof = —0.25. We show the distribution in the 
(t, r 3 ) plane so that a uniform distribution of points in the 
plot corresponds to a uniform distribution in the 3 dimen- 
sional space. The distributions are based on summing the 
results of 10 events. In the thermal explosion, a total of 362 
last scatterings took place in the 10 events, whereas in the 
implosion all particles scattered at least once in each event 
(500 scatterings in total). 

An attractive way of understanding the freeze-out pro- 
cess is to think of the field of most recent scatterings (the 
scattering field): At any given time, there will be points 
in space where each particle had its most recent scatter- 
ing. The distribution of these points in space-time is the 
scattering field Z(r, t; t'), where f and t are space-time co- 
ordinates and t' is the present time (the instant at which 
we look at the system). The entire field I as function of 
r and t will of course change with time i', but only un- 
til complete freeze-out when the scattering field becomes 
fully fixed and equal to l{r, t; 00). 

In addition to the freeze-out field itself, we can con- 
sider other variables defined on the freeze-out space-time 
coordinates. Of particular importance is the freeze-out 
flow field, since this may be related to the experimen- 
tally measured flow. The freeze-out flow field should be 
defined as a non-thermal part of the particle velocity at 
the freeze-out points. In general the definition of collec- 
tive flow is somewhat ambiguous. However, in our case 
of spherical geometry it can be defined as the mean ra- 
dial component of the particle velocities at the freeze-out 
points. In Fig. ||| we show some examples of time inte- 
grated freeze-out flow fields. Any non-zero value of these 
curves is a sign of collective flow. 
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FIG. 3. Time integrated freeze-out flow fields for 
voj = —0.25 (the three dashed curves) and vof = —0.5 (the 
three solid curves) . Within each set of curves the flow fraction 
r\ in the initial state is varied, n — 0, 0.5, 1 (bottom to top). 

The general picture is that the freeze-out flow increases 
with the ingoing flow fraction r/ and with the distance to 
the center of the explosion. Note also that even for the 
thermal explosion (j) = 0), the freeze-out flow is non-zero. 
One sees from the figure that within statistical errors the 
two families of curves reveal the scaling with the initial 
velocity. 

Equilibration. We do not wish to make any assump- 
tions about equilibrium, be it global or local, partial or 
complete. We will, however, make use of the concept of 
entropy, since this can formally be calculated also out of 
equilibrium. 

To study one-body observables, we reduce the 6 A di- 
mensional phase-space of the A particles to the one-body 
phase-space in 6 dimensions in the standard way jlO) . We 
further exploit the radial symmetry of the system by re- 
ducing the phase-space to two dimensions (|r|,|p|), the 
advantage being a gain in statistics. Then we introduce 
a finite grid in this reduced phase-space, dividing each of 
the two axes into a number of segments. But instead of 
working with a fixed grid in phase-space, we let the entire 
grid expand or contract with the majority of particles. 
Our experience shows that i? 2 o = (1M J2 n =i r 20 ) 1 / 20 is 
a good measure of the spatial extension of the system. 
This means that the physical size fii (in units of (27r?i) 3 ) 
of cell i in phase-space varies with time. We choose our 
grid such that all cells have equal phase space volume. 

One can use the phase-space distribution of particles 
obtained from microscopic simulations to characterize 
the degree of equilibrium. First of all, for given ini- 
tial conditions we perform several independent dynam- 
ical runs, and accumulate the points in different phase- 
space cells at each time step. In this way we obtain the 
occupation probability distribution {p{\, where 

number of points in box i 
Pi = total number of points in phase space ' (2) 



Then we repeat the simulations for a reference system 
where particles have sufficient time for equilibration in a 
fixed volume 47r/3i?2 , where i? 2 o is calculated from the 
actual spatial distribution of particles in the dynamical 
simulations. In this way we obtain the corresponding 
equilibrium distribution {p^ 9 }. We then construct the 
quantity 

S = ~I>l°g(PiM). ( 3 ) 

i 

which in the case of equilibrium would be the real en- 
tropy |ll| ] . Dividing by fij inside the logarithm eliminates 
the trivial trend that a larger phase space cell holds more 
points than a smaller cell. As a measure of equilibrium 
we consider the quantity 

£ = exp(S -S ref ). (4) 

£ is thus a measure of how close the actual phase-space 
distribution is to that of an equilibrized system with the 
same volume and energy. S = 1 cooresponds to full 
equilibrium. 

Since the grid expands or contracts with the particles, 
the information on flow is somewhat erased. One would 
thus expect S to reach a constant value after freeze-out. 

In Fig. ^| we show the behavior of the equilibrium mea- 
sure £ in the following three cases: (a): The particles are 
kept inside a spherical container of radius R — 4.4 fm un- 
til t = 20 fm/c, when the container walls are removed; 
(b): 100% ingoing flow {-q = 1), intended to simulate 
an implosion-explosion process; (c): The particles are 
started with 50% thermal energy and 50% outgoing flow 
(rj = 0.5), simulating an explosion from a not fully ther- 
malized state. 
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FIG. 4. The equilibration measure E calculated for the 
three cases discussed in the text. The results are based on an 
emsemble of 20 events with 50 particles. Here, r c = 0.4 fm. 

We notice the following features of £ in this figure: 
In the case (a) S ~ 1 as long as particles are kept in 
the container, reflecting the fact that the system is in 
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equilibrium. Then when the container is removed, the 
system goes out of equilibrium and £ decreases until col- 
lisions between the particles cease. This is as intuitively 
expected. Now, in the implosion-explosion case (b), we 
do not know in advance if the system will reach a state of 
equilibrium or not. Here, £ starts from a low value, since 
the initial state is far from equilibrium. As collisions tend 
to equilibrizc the system, £ increases, but does not reach 
the equilibrium £ = 1. However, additional analyses 
show that the particles obtain a nearly Maxwellian veloc- 
ity distribution from t — 6 fm/c, which is also the time 
of maximum compression. During the de-equilibration 
time, £ decreases, and the rate of decrease reflects the 
dynamics of freeze-out. It is interesting to note that the 
thermal explosion (a) and the "non-thermal explosion" 
(c) both approach the same asymptotic value £ ~ 0.55, 
as does the implosion-explosion case. Therefore we can- 
not use the asymtotic behavior of E to conclude about 
the degree of equilibrium at early stages of the system 
evolution. The asymptotic value of E is not universal, 
for instance for 100% outgoing flow E is approximately 
constant ~ 0.4. 

Conclusions. We have performed dynamical simula- 
tions within a simple gas model and analysed the re- 
sults both with regard to freeze-out, flow, and degree of 
equilibrium. The space-time distributions of last scatter- 
ings are broad in both space and time. Without using 
any assumptions about equilibrium we have defined the 
freeze-out collective flow, at least for the simple geome- 
try adopted here. A general trend is that the freeze-out 
distribution is sharper and the freeze-out flow is stronger 
when the system goes through a state of high compres- 
sion. Our calculations show that the maximum compres- 
sion coincides with the maximum degree of equilibrium, 
but even at this stage equilibration is not complete. The 
present investigation uses a very schematic interaction. 
In the future we are planning to extend the study to 
other types of interaction. 

Our general conclusion is that small systems evolving 
under fast dynamics exhibit many features which cannot 
be described adequately by equilibrium concepts. 
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